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INFINITE POTENTIAL RELAXED 
by Gabriel Allen 
Lewis Research Center 


SUMMARY 

Exact solutions have been obtained for a model representing an electron in 
a one -dimensional periodic potential acted on by a transverse uniform magnetic 
field. The periodic potential is approximated by a finite chain (20 atoms 
long) of "periodic" square wells. Within this chain, the parabolic potential 
due to the magnetic field was approximated by a parabolic square-well potential. 
The infinite potential at each end of the chain, which was used in a previous 
report to cut off the chain sharply, is relaxed and the true magnetic parabolic 
potential was used outside of the chain. Eigenvalues and wave functions are 
computed for the case of equal widths of wells and hills in the periodic part 
of the potential. A reasonably consistent interpretation of the behavior of 
the wave functions is presented with crossing-like behavior between zero field 
energy bands being a likely part of the interpretation. 


INTRODUCTION 


In a previous report (ref. l) exact quantum-mechanical solutions were ob- 
tained for a model representing an electron in a one -dimensional periodic po- 
tential acted on by a uniform transverse magnetic field. The results of an 
analysis of the eigenvalues and wave functions of this system showed that the 
model seemed to be capable of following certain magnetic breakdown effects in 
some detail. It was possible to see the increasing influence of the magnetic 
fields on the behavior of the states by following individual states from zero 
magnetic field (where the periodic band structure properties control the be- 
havior of the system) to magnetic fields larger than Blount breakdown fields. 


For the system in reference 1, the Schroedinger equation was reduced to a 
one -dimensional problem with an effective potential of the form shown in fig- 
ure 1 (called a parabolic square-well potential). In the figure, the depth of 
the periodic potential well is denoted by Vq and the magnetic step in the po- 


tential is denoted by where 


V* 

V 0 



CD 


where a is the period. 



The absence of periodicity meant that the number of matching conditions 
(and thereby the rank of the determinant used for the determinantal compati- 
bility condition) increased with the length of the chain. The absolute cutoff 
condition, which may be written as 

V |(1T + l)a - wj ■ co (2) 

was imposed to keep the number of computations down to a manageable level. In 
equation (2), N is the number of atoms in one-half of the chain, and 2w is 
the width of the periodic potential well. 

The main effect of this condition appeared to be that wave functions be- 
longing to eigenvalues near Vq + were forced to zero near the cutoff 

region at the end of the chain. The relaxing of the absolute cutoff condition 
should therefore improve the model. To this end, a program for determining 
eigenvalues and wave functions for an effective potential of the form shown in 
figure 2 was undertaken. The program was completed and the results have been 
used to try to interpret those effects in both the eigenvalues and wave func- 
tions arising from the cutoff condition. 

Significant changes in the results were not expected, and the main fea- 
tures of the preceding model were found again. However, the effects of the 
possible interchange of properties between different bands were not considered 
in reference 1, and the likelihood of such events permitted a more detailed 
interpretation of the results to be made than in the preceding report. 


SYMBOLS 

A m constant multiplying wave function in pure magnetic region 

a period 

c velocity of light 

functional form of determinantal compatibility condition 
D v (|) v th order parabolic cylinder function 

E parametrized quantity used like energy, see eq. (6) 

e charge of electron 

G 

g defined by eqs. (ll) and (l2) 

H magnetic field strength 
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H 


magnetic field 


Eg Blount breakdown value of magnetic field 

H z magnitude of z -component of magnetic field 

ft Dirac h, Planck* s constant/2it 

k z component of wave vector 

m mass of electron 

N number of atoms in one -half of chain 

n integer defining level of state in hand 

r cyclotron radius 

B-b; defined by eqs. (All) and (A12) 

B^sJ defined by eq. (A13) 

S^.S defined by eq. (A14) 

n n 

Vq well depth of periodic potential 

Vq magnetic step in potential 

V(x*) total potential 

w half -width of periodic potential well 

x,,y,z coordinates 

e actual eigenvalue of system 

A(x T ) wave function 

v defined by eq. (l5) 

£ defined by eq. (l6) 

cd c cyclotron frequency 
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Subscripts : 

M pure magnetic 

N last period in periodic potential region 

n n^* 1 period 
Superscripts : 
h hill 
w well 
R reduced 

0 sign of appropriate g 


MODEL DESCRIPTION 


Most of the characteristics of the model are the same as in reference 1. 
In particular, in the region 


[(N + l)a - wj < x' < |(N + l)a - wj 


the same parabolic square-well potential applies as in reference 1. (Compare 
this region in figs. 1 and 2.) However, in reference 1, at 


x r = (N + l)a - w 

the cutoff condition given by equation (2) is applied. This condition may be 
described as an absolute cutoff in the sense that it has the effect of keeping 
the electron entirely inside the "box" so that for 

| x f | > (N + l)a - w 

the wave function A(x r ) is zero. 

In this report, the absolute cutoff condition is partially relaxed (by 
keeping the "pure" magnetic potential) to allow the system to penetrate into 
the region |x T | > (N + l)a - w. However, the periodic part of the potential 
is still cut off in order to keep the number of matching conditions finite. In 
this region, the effective potential will subsequently be called pure magnetic. 

Relaxing the absolute cutoff condition makes for an improved model in the 
obvious sense of releasing the electron from its "box. 11 However, the boxing-in 
aspect had the advantage of providing a natural starting point for the individ- 
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Figure 1. - Form of potential given by parabolic square well with absolute cutoff. 



Figure 2. - Form of potential due to parabolic square well without absolute 
cutoff, w = a/4. 

— > 

ual states at H = 0. The identity of these zero-field states was preserved 
and the individual states could be followed as h increased to Blount break- 
down values . 


In the present model, the absolute cutoff condition is replaced by the 
pure magnetic potential 


V(x') 


l (r“i x ’ 2 J 


in the region |x'| > (W + l)a - w, where the cyclotron frequency is 


When figure 2 is referred to, it can be seen that if H = 0 there will be 
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no bound states at all; however, any finite H, no matter how small, will bind 
the system. Thus, there can be no easy transition from small fields to zero 
field. 

There is no good reason to give up this advantage of the cutoff model. For 
the lower lying states, the absolute cutoff condition has very little effect on 
the system since for the fields studied in reference 1 the finite potential bar- 
rier was already very large at x r < (1ST + l)a - w. Thus, the main advantage 
gained by relaxing the absolute cutoff was in examining the behavior of higher 
states (at comparatively small magnetic fields) whose amplitudes would still be 
appreciable in this vicinity if not forced to zero by the cutoff. (At high mag- 
netic fields, the cyclotron radius becomes so small that the wave function has 
a very small amplitude near the periodic cutoff except for states at rather high 
levels.) Now it is desirable to follow the behavior of states which flank the 
first energy gap (band gap) at zero field up through higher fields. At rela- 
tively low fields, those states would come under the heading of low-lying states 
as described previously. There is thus no reason to relax the absolute cutoff 
condition at small fields and the eigenvalues for zero field computed using this 
condition will be retained in this model. 

A proper way to effect a smooth transition would be to start with a fully 
periodic potential (a "chain" of infinite length) at zero field and then intro- 
duce a periodic cutoff potential of monotonically decreasing length as H in- 
creases. Computational difficulties make such a procedure impractical. 

To summarize then, the model^used herein is the parabolic square well with 
periodic cutoff of figure 2 for H > 0 and the periodic potential in a box of 
figure 1 at H * 0. 


SOLUTION OF WAVE EQUATIONS 

After the standard reduction to one dimension, the wave equation for the 
physical system becomes 


d^A 


dx' 



where V(x*) is the potential shown in figure 2 and 



( 5 ) 


( 6 ) 


where e is the actual eigenvalue of the system and k z is that component of 
the wave vector k which is parallel to the magnetic field H. These equa- 
tions and definitions are identical with equations (23) and (2l) respectively in 
reference 1, except for the fact that in this report the cutoff condition (2) is 
replaced by the condition 
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V(x,) = ^(^ mCD c x ' 2 ) l X 'l 


> (N + l)a - w 


( 7 ) 


The solution to equation (5) in the regions of constant potential contained 
in the interval 0 < x' < Na + w are of the same form as in reference 1 and are 
summarized below for the even solutions: 


, w> t \ ^w „w , 
AqCx 1 ) * B Q cos G q x' 


( 8 ) 


where -w < x* < w, 

r 


— sinh G^ - (na - v)J + B n cosh G n - (na - w)J g n < 0 
£(*') “ A n [ x ' " ^ na " w ^] + B n 


V 

y 

3 n 


w 


sin G^ jx' - (na - w)J + cos G^jx' - (na - w) 


g W > 0 

B n 


n 


( 9 ) 


where na - w < x' < na + w, and 


r 

sinh G^ Jx* - (na + w)J + B^ jx’ - (na + w)J 
G. 


h 


h 


n 


A n (x’) = < 


jx‘ - (na + w)J + B^ 
h 

sin G^ jx* - (na + w)j + B^ cos G^ |x* - (na + w)J 


h 

g n < 0 


«n “ 0 


h 

^n 


> 0 


n 




where na + w < x* < (n + l)a - w. In equations (8) to (10), 




2m 


i,2, 2 N 

k z 2. • 

"ST “ n V 


( 10 ) 


( 11 ) 
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( 12 ) 


g 


h 

n 



fi 2 k| 

2m 



g W “ 
s n 


2m 




2 V 0 


G = 



( 13 ) 


Note that in this report equation (10) also applies in the region of the 
last hill, Na + w < x’ < (N + l)a - w, unlike the case in reference 1 which re- 
quires a special form for this region. 


For x* > (N + l)a - w, the harmonic oscillator potential i mu& x ' 2 
governs the behavior of the system and the solution becomes 


a m ( X ') - AmD v U) 


( 14 ) 


where D y (£) is the order parabolic cylinder function (ref. 2), is a 

constant, and 


v 


€ 1 

hu> “ 2 
c 





X * 


( 15 ) 

( 16 ) 


The subscript M denotes the fact that in this region the potential is pure 
magnetic. 


The modifications that this requires in the matching conditions and their 
use in computing eigenvalues are shown in appendix A. The use of Weber func- 
tions of high order introduced some special difficulties in the computations, 
which are discussed in appendix B. As in reference 1, the eigenvalues actually 

2 2 

computed are e - h k z /2m rather than e. 


WAVE FUNCTIONS 

Wave functions were obtained in the usual way by finding the coefficients 
in each region of constant potential in terms of an arbitrary coefficient which 

was the coefficient in the well, B^. The procedure was essentially the 
same as described in reference 1 and need not be repeated here. 

It should be noted that the wave functions were not normalized in the 
usual fashion. The difficulty lay in the fact that a proper normalization 
would require the evaluation of integrals of the form 

f DyU) 

e M 
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where 


°M 


/SmaiA 1 / 2 r n 

■(— ) [ (M+1)a - w ] 


There is no easy way to evaluate such integrals. Therefore, an incomplete 
normalization was used. The procedure consisted merely of neglecting the inter- 
val in the aforementioned integral and using the condition 




(N+l)a-w 


|A(x')| 2 dx* = i 


(16) 


On the one hand, the use of this condition means that quantitative compar- 
isons of the amplitudes of wave functions belonging to different eigenvalues 
are not valid. On the other hand, the kind of comparison that is suitable for 
the particular analysis desired can be made quite rigorously. For a given 
state, the amplitude is examined as a function of x* to determine in which 
part of the "crystal” the system has a comparatively larger probability of 
being found. Then, a similar examination is made for a different state in 
order to determine whether the two states have large and small amplitudes in 
the same regions. In this connection, it may be mentioned that for most states 
the neglected contribution to the normalization is quite small so that even the 
approximation to the absolute magnitudes of the amplitudes is good. 


RESULTS AM) DISCUSSION 
Eigenvalues 

The eigenvalues for various well depths are shown as a function of mag- 
netic field strength H in figure 3. Since magnetic field strength is a com- 
mon parameter, it was chosen as the ordinate in preference to V Q . However, 

V* is also shown on the vertical scale to make it easier to compare these re- 
0 

suits with those in reference 1. The discontinuities in the slopes in the fig- 
ures are due to the procedure used. For each V computations were carried 

out for a few values of H and these eigenvalues were then connected by 
straight lines. It is convenient to let these lines represent the same "state" 
at different magnetic field strengths. As explained in the last section, the 
eigenvalues for H z = 0 are those taken from the box model of reference 1. 

The Blount breakdown lines are shown on the figures as an aid in the analysis 
of the results. General quantum mechanical considerations ensure that alternate 
even and odd solutions occur (ref. 3). In a few cases, D^ oscillated near but 

not through zero for a series of trials for even solutions in an energy inter- 
val in which several odd solutions were found. In these instances, circles 
were placed between odd eigenvalues in the general vicinity of where the even 
eigenvalues must lie. 
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A fact should be noted here that was perhaps not made sufficiently clear 
in reference 1 to which it also applies. In computing the eigenvalues , an at- 
tempt would be made to find a value of e which would make (see appen- 

dix A) smaller than some preassigned value. This value of Djj which was chosen 
by trial and error was small enough to enable the boundary conditions to be sat- 
isfied by the wave functions. Although only those eigenvalues which satisfied 
this criterion were used in computing the wave functions, all of the values of 
€ at which % changed signs are presented in figure 3 both here and in ref- 
erence 1. The identification of the band to which a given eigenstate belongs, 
together with its location in the band, is important even if the eigenvalue it- 
self cannot be used in computing a satisfactory wave function. 

Wave Functions 

Before discussing the results of the calculations in this report, it might 
be of some interest to review briefly the results in reference 1. The general 

(where Hg is the 

Blount breakdown value) the wave functions were greatly influenced by the peri- 
odic part of the potential. Thus the relative probability that the system 
would be found near compared to the probability that it would be found far from 
the cutoff was quite similar for states within the same band but changed dras- 
tically in going from a state on one side of the energy gap to a state just the 
other side of this gap. This, change in behavior of the eigenstates on crossing 
an energy gap was always great even when the difference in the energies of the 
states flanking the gap at a fixed value of H was very small (see figs. 3(c) 
and 5 in ref. l) . As the magnetic field increased, however, the effect on 
the wave functions of crossing the energy gap became less pronounced -until, at 
H ^ Bg, no discernible effect remained (see fig. 10, ref. l) . 

It may be argued that the use of negative amplitudes was unnecessary. The 
reason for plotting the actual wave function rather than | A | or | A | ^ was that 
it was easier to check visually the continuity of the slope across boundaries 
separating regions of different potential from A(x f ) than from either of the 
other choices. 

When describing the strong effect of the periodic part of the potential on 
the behavior of these states, there is no intent to imply that any departure 
from the general behavior of electrons under the combined influence of periodic 
potentials and magnetic fields is being shown by this model. The magnetic 
field is the principal influence on the electron (ref. 4). All that the model 
is demonstrating is the difference between the effect of a magnetic field on 
free electrons and on electrons in a periodic potential. Electrons in states 
near the bottom of a band act very much like free electrons. Wave functions 
for such states are peaked near their cyclotron radius 

r = jfrk _ 11,5 V n + 2 

where is in kilotesla and n is the number identifying the level of the 


behavior of the wave functions was such that for H « Eg 
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state in the band (n = 1 for the ground state) . For the model used here, these 
peaks are near the center of the effective potential. 

As the states approach the top of the band, their behavior becomes less 
free and is more difficult to describe. However, at the bottom of the next 
band, a "free electron"-like behavior is to be expected again with peaks rela- 
tively near the center of the effective potential. For very large values of 
H, the effect of the periodic potential becomes less pronounced and the behavior 
of states near the top of a band will not be very different from those near the 
bottom of the successive band. 

The preceding description omits any discussion of the possibility that the 
magnetic field can cause a crossing between states originally in different 
bands. Such an omission appears quite acceptable at first sight. If a given 
solution to a one -dimensional Schroedinger equation is a continuous function of 
a parameter, then the number of nodes of the solution will not be changed by a 
change in the parameter. Now the magnetic field H z is a parameter of the 
wave function here, so a wave function that has, for example, 20 nodes when 
H = 0 will have 20 nodes when H = 5 kilotesla. Since the state represented 
by this function was in the first band at H = 0 (the state at the top of the 
first band has 21 nodes), the state which has 20 nodes when H = 5 kilotesla 
still has to be regarded as the same state (insofar as its location in the band 
is concerned) . 

This fact does not prevent the properties of the wave function at large 
H z from being very much like those from another band at small H z . A qualita- 
tive description of the situation can be made with the aid of figure 4(a), 
which is a schematic representation of the variation of the two distinct states 
1 and 2 with an arbitrary parameter p. At a certain value p c of the param- 
eter, there is a near crossing of states 1 and 2. Under such circumstances, 
the wave functions for either of the two states may acquire a large admixture 
of properties of the other state. If such a mixing of properties occurs, then 
at some larger value of the parameter, for example p Q , the wave function for 
state^ 1 might have properties much like those which would be expected from an 
extrapolation of properties of state 2 when p « p c and vice versa. If 
states 1 and 2 are in different bands for p « p c , the aforementioned pattern 
of events will be called an "interband transfer." 

Figure 4(b) shows a possible two-stage interband transfer between a pair 
of states originally in the first band and a second pair originally in the sec- 
ond band. 

The results of the present computations are shown in figures 5 and 6 for a 
Vq of 5 and 10 electron volts, respectively. In examining these figures, it 
should be kept in mind that the normalization was not complete so that quanti- 
tative comparisons between the amplitudes of different states at the same posi- 
tion may not be reliable. On the other hand, for a given state, the ratio of 
the amplitude near the periodic cutoff to the amplitude near the center will be 
independent of the absolute magnitude of the amplitude. A comparison of such 
ratios between two different states will, therefore, be reliable. 
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(a) Schematic possible transfer of properties between two states having a near-crossing. 



(b) Possible two- stage interband transfer between two states at top of first band and two states 
at bottom of second band. Transfer occurs between states ly and 2g at (a), between ly and 
ly at (b), between ly and 2g again at (d), and between 2g and 2g at (c). For large values 
of magnetic field H, ly and ly states would have properties like states 2g and 2g, 
respectively, while 2g and 2g would have properties like ly and ly, respectively. (State at 
bottom of first band, lg; state at top of first band, ly; state just below top of first band, If; 
state just below state ly, If; state at bottom of second band, 2g; state just above bottom of 
secondhand, 2g; state just above state 2g, 2g*. ) 

Figure 4. - Schematic representation of interband transfer. 
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In the remainder of the discussion the expression " state in the 

hand" will mean "wave function corresponding to the eigenvalue in 

the hand. " Both the hand and the position in the hand refer to zero 

field. 

Figure 5(a) may he selected for a detailed discussion. As can he seen 
from figure 3(a ) , a field strength of 3 kilotesla (« H-g) for V 0 = 5 electron 

volts. Thus, the influence of the periodic part of the potential should he 

discernible in figure 5(a). Examination of this figure shows that interband 
transfer is evident. The tenth odd state in the first band behaves like a 
state nearer the bottom of a band than the top. Such would be the case if the 
characteristics of states near the bottom of the second band crossed with the 
characteristics of some of the states near the top of" the first band. The next 
two higher states (which, by the way, correspond to the top of the first band 
and bottom of the second band, respectively) act like states near the top of a 

band, which is in agreement with the description of transfer behavior. The 

next two states are both second band states. They do indeed behave like second 
band states and furthermore seem to behave like states not quite as near to the 
bottom of a band as the tenth odd state. The entire set of four states shown 
thus seems to be in complete conformity with the interpretation of interband 
transfer described in the preceding paragraphs. 

Next figure 5(b) may be examined. The interpretation here is also in 
agreement with the overall picture described, but the results are less striking. 
The tenth odd state decidedly shows the behavior expected from a state near the 
top of a band so the crossing, in evidence for this state at 3 kilotesla, has 
disappeared at 4 kilotesla. The behavior of the other three states shown on 
the figure does not admit of a clear interpretation in terms of band position 
since no definite preference is shown for any particular part of the atomic 
chain. A comparison with figure 10(a) in reference 1 can be made to attempt to 
ascertain the effect of the absolute cutoff. That figure showed the behavior 
of wave functions in a potential having Vq = 5 electron volts and H = 3.8 
kilotesla, which is close to the conditions holding for figure 5(b) in this re- 
port. 


It may be noted that the cutoff states seem to be showing transfer behav- 
ior since the states near the top of the first band at zero field exhibit the 
behavior of states near the bottom of a band, whereas the state which repre- 
sented the bottom of the second band at zero field behaves like a state near 
the top of a band. The absolute cutoff thus seemed to effect an interband 
transfer that is not present when the cutoff condition is relaxed. 

Figures 5(c) and (d) can now be examined together and compared with fig- 
ures 10(b), (c), and (d) in reference 1. The potential for the state in fig- 
ure 5(c) of this report is almost exactly the same as for the states in figure 
10(b), whereas the magnetic field in 5(d) is about midway between those in fig- 
ures 10(c) and (d) of reference 1. The suggested comparison shows that the be- 
havior of the states is changed very little by relaxing of the cutoff condition. 
As asserted in the description of the model, at high magnetic fields, the ef- 
fect of the cutoff should not be very strong, and this is precisely what the 
results confirm. In particular, the phenomenon of Blount breakdown and the 
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(d) Magnetic field strength, 6.0 kilotesla (near Blount breakdown). 
Figure 5. - Concluded. 
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disappearance of the effect of the periodic' potential is retained upon the re- 
laxation of the cutoff condition. 

The discussion of the Vq = 10 case may he prefaced by some pertinent re- 
marks. First, improved computations enabled a more complete set of wave func- 
tions to be presented here than was available in reference 1. However, some 
eigenvalues still did not have sufficiently small values to permit their 

use in the computation of wave functions. The fact that the set of wave func- 
tions for Vq = 10 was relatively incomplete in reference 1 makes it difficult 
to compare the wave functions with and without the cutoff for this case. What 
comparison is possible, however, is in agreement with the general ideas pre- 
sented herein. Secondly, the larger zero field energy gap would tend to re- 
quire higher magnetic fields to get interband transfer and other wave function 
behavior similar to the Vq = 5 behavior. 

Figure 6(a) (H = 3.0 kilotesla) exhibits the effect of the cutoff very 
strikingly. A comparison of this figure with figure 5 of reference 1 (which 
shows wave functions at H = 3.1 kilotesla) shows that the cutoff indeed forced 
the wave functions to zero near the end of the chain. This effect was so strong 
that it masked other effects (such as band effects) which might otherwise have 
been manifest. Also, the eigenvalues themselves were pushed noticeably higher 
(the states flanking the energy gap were displaced upwards.) 

As shown in figure 6(b), when H = 4.0 kilotesla, the wave functions behave 
in a manner qualitatively similar to those in figure 5(a). The 13.91-electron- 
volt state is clearly in the second band, the next two are back in the upper 
part of the first band, and finally the second band states continue to fill up. 

Figures 6(c) to 6(f) show the wave functions as H goes from 5.0 to 
8.0 kilotesla. The same general behavior is shown by the eigenstates as for 
the shutdown V 0 = 5 .0-electron-volt potential wells. Interband transfers 
show up rather well at the lower magnetic fields; then, as H increases, 
the evidence of the band structure gradually disappears until at 7.0 kilotesla 
no remnant of this structure remains. It should be noted that the effects of 
the periodic part of the potential became small at values of H considerably 
less than H-g (see fig. 3(b)). This phenomenon might not be quite the same as 
Blount breakdown; it may rather be a precursor of it. 


SUMMARY OF RESULTS 

A one -dimensional model has been examined to try to follow some details of 
the magnetic breakdown process. The model is that of an electron in a uniform 
magnetic field H z and a one-dimensional transverse chain of 20 periodic 
square wells. This model is an improvement over a previous one which cut off 
the chain by means of an infinite potential at each end. 

Exact solutions were obtained for both the eigenvalues and wave functions. 
Examination of the wave functions showed that, for comparatively small magnetic 
fields, at a fixed magnetic field and periodic well depth Vq the behavior of 
the eigenstates admits a rather consistent interpretation in terms of an inter- 
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band transfer of characteristics between bands. As the magnetic field is in- 
creased, however, the band effects become less pronounced. For Vq = 5 elec- 
tron volts, the value of the magnetic field at which such band structures ef- 
fects become unimportant is the Blount breakdown value H-g. For the deep wells 
(Vq = 10 eV) , the magnetic field strengths are still quite a bit lower than the 
Hg values, which indicates that the system may be in some condition which is a 
sort of precursor of an actual magnetic breakdown state. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, October 22, 1965. 
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APPENDIX A 




MATCHING CONDITIONS AND DERIVATION OF EIGENVALUES 

The matching conditions up to the N^* 1 hill are the same as those in refer- 
ence 1 and are found in appendix A therein. In the N't* 1 hill, however, 

— sinh G^J (x 1 - (Na + wj] 
x' - (Na + w) 
sin g n [ x ’ ~ ( Wa + w 0 


G 


N 


(x* ) - A h J 

N N ' 




G: 


N 


+ < 


cosh G^ Jx' - (Na + wjj 
1 

cos G* 1 fx' - (Na + wj] 

TvT 1 — -J 




4 <o 

} ■ 4- 0 Cal) 

g h > 0 
N 


and 


^cosh G^ Jjx 1 - (Na + wjj 


^ = A* 1 ^ 

dx’ V ‘ n ’ 




cos G^ |x' - (Na + wjj 


+ B h G h j 
N N ' 


-\ 


sinh G* 1 |x' - (Na + wjj 


-sin G^ Qc* - (Na + wjj 


g h < 0 

N 


> (A2) 


g h > 0 

S N 




In the pure magnetic region, x T > (N + l)a - w, 

and from reference 2, 

afr Am(x') = [§ D v (i) - D V+1 U)] 


(A3) 


(A4) 
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The matching conditions then become 


+ b M - va 


❖e*!} + e h c e b h = VA 


C L A M = D V 


C L% “ 


2m % Rl 


r + l)a - w \ = D v (i M ) 

(A7) 

D v^M^ " •’-’v+I^mU 

(A8) 

(N + l)a - wj 

(A9) 


and the other quantities appear in appendix A of reference 1. 

With these changes, the determinantal compatibility condition is obtained 
by setting equal to zero the following determinant which corresponds to 

equation (A3l) in reference 1 (note that this is specifically for the case of 
even functions) : 


M "N “H “N "N-l hj-1 


-Vm °A VS 0 0 

^ ^ 0 0 


0 - 1 c r a w Cr^ 0 


0 VS C R^ 


0 0 0 


0 0 0 


0 0 
0 0 


- 1 C R A N-1 


0 C X-1 C R B tl 


0 0 0 0 0 0 


0 0-1 


0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
0 0 


0 0 
0 0 


Vo Vo 


-1 0 -G" sin wG" 

0 0 


fin ^nui ■■■mi mi iiiii ii ii ii mi 1 111111111 11 1 1 11 nun 11 


1 11 hi 11 1 11 11 1 111 mini min ii 111 iniB 1 in 11 ■■ 11 mill 1 1 111 1111111 111111111 1111111 111 11 1 1 1 1 m m in iiihiii i iiiiiiii iiiiiiiiii ii 1111 hi 



The procedure for evaluating is the same as was used in 

but will be sketched here for completeness. First the quantities 
are defined as follows: 


For even solutions: 


For odd solutions: 


riW ~w 

Sq = cos wGq 




S£ = -GY sin wG v 


0 


0 


rtW 1 . «W 

= — sin wG_ 
0 w 0 

G 0 


Sq = cos wGq 


Then equations (A13) are used alternately with equations (A14) up 
eluding n = R: 


where 


S h = CL A h s w + C T1 B h S w 
n R n n R n n 

-h r h— w f h w f 

S = C d A S +C d BS 
n R n n R n n I 


S “ - + °E® 


! W S h A 
n n-ll 


S w = C A"S 


iWTTh. 


n 


R n n-1 


+ CLB W S 


,w 0 h 


R n n-lj 


C+A w 
R n 


- < 


— sinh 2wG w 
w n 

G n 


2w 


— sin 2wG w 
n 


reference 1 

s o G o 


(All) 


(A12) 


to and in- 


(A13) 


(A14) 


(A15) 
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Finally, 


cosh 2wG. 


0 0 
n^viy ~~ p * — 

TTn U R " 


cos 2wG. 


,w 

-n 


0 

C + B w = < 
R N ^ 


r 


G W sinh 2wG W 
n n 


~G sin 2wG 


n 


0 

C*A h = <( 
R n ^ 


1 . , 

— sinh 2hG 

G n 

2h 


sin 2hG h 

o h 

n 


n 


0 0 

C + B h = c’ + A h =* 1 
R n R n ' 


cosh 2hG; 


h 

~n 


cos 2hG 


n 


r G^ sinh 2hG^ 


0 

C ' + B h = < 
R n 


0 

-G h sin 2hG h 
n n 


D » - -S^l/Sd + S S C L*M - 0 


(A16) 


(A17) 


(A18) 


(A19) 


( A20) 


(A21) 


These quantities are used in the computation of D^. When is deemed 

to be sufficiently close to zero, the value of E which was used in that com- 
putation is considered to be an eigenvalue for the given combination of Vq 

and H z • 
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APPENDIX B 


COMPUTATION OF WEBER FUNCTIONS 

Parabolic cylinder functions (or Weber functions) oscillate with amplitudes 

proportional to V* + “•). From equation (14) it can be seen that v can be 

very large in these computations and, in fact, became large enough so that the 
numbers in the computations would sometimes fall outside the range of the Lewis 
7094. For this reason, a "reduced" Weber function with an amplitude not exceed- 
ing one unit was used in this computation (ref. 6). A description of the use 
of these functions follows. 


The superscript R denotes reduced. Then the determinantal compatibility 
condition = 0 can be written as 

d n = + |) D N = 0 (B1) 

It may be noted that for the values of E in this report, r(v + \)> "will always 
be greater than zero so that and vanish at the same values of E. 

Here 

4 = + < E2 > 


where 


c'a r 

11 


c'a r 

°l a m 




(B3) 


V 2mco | ”] 


(B4) 


Practical procedures for computing values of high-order Weber functions do 
not seem to be readily available especially at large arguments (ref. 7), so a 
procedure for obtaining such values on the Lewis 7094 was developed. 
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